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Abstract. Consider a deterministic self-adjoint matrix X n with spectral measure con- 
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ues converging to the edges of the limiting measure. We perturb this matrix by adding 
a random finite rank matrix with delocalised eigenvectors and study the extreme eigen- 
^-j- values of the deformed model. We give necessary conditions on the deterministic matrix 

X n so that the eigenvalues converging out of the bulk exhibit Gaussian fluctuations, 
whereas the eigenvalues sticking to the edges are very close to the eigenvalues of the 
non-perturbed model and fluctuate in the same scale. 

We generalize these results to the case when X n is random and get similar behavior 
when we deform some classical models such as Wigner or Wishart matrices with rather 
general entries or the so-called matrix models. 
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1. Introduction 



Most of the spectrum of a large matrix is not much altered if one adds a finite rank 
perturbation to the matrix, simply because of Weyl's interlacement properties of the 
eigenvalues. But the extreme eigenvalues, depending on the strength of the perturbation, 
can either stick to the extreme eigenvalues of the non-perturbed matrix or deviate to some 
larger values. This phenomenon was made precise in [9], where a sharp phase transition, 
known as the BBP transition [3U E7J |3U EH], was exhibited for finite rank perturbations 
of a complex Gaussian Wishart matrix. In this case, it was shown that if the strength 
of the perturbation is above a threshold, the largest eigenvalue of the perturbed matrix 
deviates away from the bulk and has then Gaussian fluctuations, otherwise it sticks to the 
bulk and fluctuates according to the Tracy- Widom law. The fluctuations of the extreme 
eigenvalues which deviate from the bulk were studied as well when the non-perturbed 
matrix is a Wishart (or Wigner) matrix with non-Gaussian entries; they were shown to 
be Gaussian if the perturbation is chosen randomly with i.i.d. entries in [7], or with com- 
pletely delocalised eigenvectors [18j fl9~] , whereas in [12], a non-Gaussian behaviour was 
exhibited when the perturbation has localised eigenvectors. The influence of the localisa- 
tion of the eigenvectors of the perturbation was studied more precisely in |13j . 

In this paper, we also focus on the behaviour of the extreme eigenvalues of a finite rank 
perturbation of a large matrix, this time in the framework where the large matrix is deter- 
ministic whereas the perturbation has delocalised random eigenvectors. We show that the 
eigenvalues which deviate away from the bulk have Gaussian fluctuations, whereas those 
which stick to the bulk are extremely close to the extreme eigenvalues of the non-perturbed 
matrix. In a one-dimensional perturbation situation, we can as well study the fluctuations 
of the next eigenvalues, for instance showing that if the first eigenvalue deviates from the 
bulk, the second eigenvalue will stick to the first eigenvalue of the non-perturbed ma- 
trix, whereas if the first eigenvalue sticks to the bulk, the second eigenvalue will be very 
close to the second eigenvalue of the non-perturbed matrix. Hence, for a one dimensional 
perturbation, the eigenvalues which stick to the bulk will fluctuate as the eigenvalues of 
the non-perturbed matrix. We can also extend these results beyond the case when the 
non-perturbed matrix is deterministic. In particular, if the non-perturbed matrix is a 
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Wishart (or Wigner) matrix with rather general entries, or a matrix model, we can use 
the universality of the fluctuations of the extreme eigenvalues of these random matrices to 
show that the pth extreme eigenvalue which sticks to the bulk fluctuates according to the 
pth dimensional Tracy- Widom law. This proves the universality of the BBP transition at 
the fluctuation level, provided the perturbation is delocalised and random. 
The reader should notice however that we do not deal with the asymptotics of eigenvalues 
corresponding to critical deformations. This probably requires a case-by-case analysis and 
may depend on the model under consideration. 

Let us now describe more precisely the models we will be dealing with. We consider 
a deterministic self-adjoint matrix X n with eigenvalues A" < • • • < A™ satisfying the 
following hypothesis. 

Hypothesis 1.1. The spectral measure /i n := n' 1 Ym=i of X n converges towards a 
deterministic probability measure fix with compact support. Moreover, the smallest and 
largest eigenvalues of X n converge respectively to a and b, the lower and upper bounds of 
the support of fix ■ 

We study the eigenvalues A" < ■ ■ • < A™ of a perturbation X n := X n + R n obtained 
from X n by adding a finite rank matrix R n = YH=i @iU™u™* . We shall assume r and the 
^j's to be deterministic and independent of n, but the column vectors (u™)i<j< r chosen 
randomly as follows. Let v be a probability measure on R or C satisfying 

Assumption 1.2. The probability measure v satisfies a logarithmic Sobolev inequality, is 
centred and has variance one. If v is not supported on we assume moreover that its 
real part and its imaginary part are independent and identically distributed. 

We consider now a random vector v n = t^(^i, • • • ,x n ) T with (xj)i<j< n i.i.d. real or 
complex random variables with law v. Then 

(1) Either the u"'s (i — 1, . . . , r) are independent copies of v n 

(2) Or (uf)i<i< r are obtained by the Gram-Schmidt orthonormalisation of r indepen- 
dent copies of a vector v n . 

We shall refer to the model (1) as the i.i.d. model and to the model (2) as the orthonor- 
malised model. 



Before giving a rough statement of our results, let us make a few remarks. 
We first recall that a probability measure v is said to satisfy a logarithmic Sobolev in- 
equality with constant c if, for any differentiable funtion / in £ 2 (V), 



/ /2log 7^^- 2c / l/,|2c ^ 



SPdi 

It is well known that a logarithmic Sobolev inequality implies sub-gaussian tails and con- 
centration estimates. The concentration properties of the measure v that will be useful 
in the proofs are detailed in Section 62 of the Appendix. 

In the orthonormalised model, if v is the standard real (resp. complex) Gaussian law, 
( u ?)i<i<r follows the uniform law on the set of orthogonal random vectors on the unit 
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sphere of M. n (resp. C n ) and by invariance by conjugation, the model coincides with the 
one studied in [TO]. 

For a general probability measure v, the r i.i.d. random vectors obtained are not nec- 
essarily linearly independent almost surely, so that the orthonormal vectors described in 
(2) are not always almost surely well defined. However, as the dimension goes to infinity, 



they are well defined with overwhelming probability when v satisfies Assumption L2 
This means the following: we shall say that a sequence of events (E n ) n >i occurs with 
overwhelming probability^] if there exist two constants C, 77 > independent of n such that 



F(E n ) > l-Ce~ n \ 

Consequently, in the sequel, we shall restrict ourselves to the event when the model (2) 
is well defined without mentioning it explicitly. 

In this work, we study the asymptotics of the eigenvalues of X n outside the spectrum 
oiX n . 

It has already been observed in similar situations, see [9J , that these eigenvalues converge 
to the boundary of the support of X n if the #j's are small enough, whereas for sufficiently 
large values of the ^'s, they stay away from the bulk of X n . More precisely, if we let G^ x 
be the Cauchy-Stieltjes transform of fix, defined, for z < a or z > b, by the formula 

G i*x( z ) = I — ^— d/xxOr), 

J Z X 

then the eigenvalues of X n outside the bulk converge to the solutions of G flx (z) = 9~ if 
they exist. 

Indeed, if we let 

o--= r ^r>o, e_-.= ^T<0 

\uxv zib G^ x {z) hnvj-a G Mx (z) 

and 

'G- x (l/9) if0e(-oo,0)U(0,+oo), 
a iiOe [6,0), 

b if 9 e (0,9), 

then we have the following theorem. Let tq e {0, . . . , r} be such that 

Ox < ■ ■ ■ < 9 ro < < 9 ro+1 <■■■ <9 r . 



Theorem 1.3. Assume that Hypothesis \1.1\ and Assumption are satisfied. For all 
i £ {1, . . . , r } ; we have 

and for all i G {r + 1, . . . , r}, 

Tn a - s \ 
A n-r+i ► Pdi- 

Moreover, for all i > r (resp. for all i > r — r ) independent of n, 

A™ ^ a (resp. \ n n ^ ^ b). 



1 Note that this is a bit different from what is called overwhelming probability by Tao and Vu but will 
be sufficient for our purpose. 
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The uniform case was proved in [TP] Theorem 2.1] and we will follow a similar strategy 
to prove Theorem 1.3 under our assumptions in Section [2] 



The main object of this paper is to study the fluct uations of the ex tre me ei genvalues 
of X n . Precise statements will be given in Theorems 3J2 3.4 , 4.3[ 4.4 and 4J3 For any 



x such that x < a or x > b, we denote by I x the set of indices i such that pg. = x. The 
results roughly state as follows. 

Theorem 1.4. Under additional hypotheses, 

(1) Let a\ < ■ ■ ■ < a q be the different values of the 9i 's such that pq. ^ {a, b} and 
denote, for each j, kj = \I Pa \ and qo the largest index so that a qo < 0. Then, the 
law of the random vector 

v ' i<i<ao v ' 9o+i<i<(? 

converges to the law of the eigenvalues of {c aj Mj)i<j< q with the Mj's being in- 
dependent matrices following the law of a kj x kj matrix from the GUE or the 
GOE, depending whether v is supported on the complex plane or the real line. The 
constant c a . is explicitly defined in Equation (|6]) . 



(2) If none of the 9{ 's are critical (i.e. equal to 9_ or 9), with overwhelming probability, 
the extreme eigenvalues converging to a orb are at distance at most n~ l+t of the 
extreme eigenvalues of X n for some e > 0. 

(3) If r = 1 and 9\ = 9 > 0, we have the following more precise picture about the 
extreme eigenvalues: 



• If Pq > b, y^A™ — pe) converges towards a Gaussian variable, whereas 
n 1_€ (A^_j — \n-i+i) vanishes in probability as n goes to infinity for any fixed 
i > 1 and some e > 0. 

• If Po = b and 9^9, n 1 ~ t (\™_ i — A„,_j) vanishes in probability as n goes to 
infinity for any fixed i > 1 and some e > 0. 

• For any fixed j > 1, ra 1_e (A™ — Aj) vanishes in probability as n goes to infinity 
for some e > 0. 

These different behaviours are illustrated in Figure [T] below. 

The first part of this theorem will be proved in Section [3j whereas Section [4] will be 
devoted to the study of the eigenvalues sticking to the bulk, i.e. to the proof of the second 
and third parts of the theorem. 

Moreover, our results can be easily generalised to non-deterministic self-adjoint matrices 
X n that satisfy our hypotheses with probability tending to one. This will allow us to study 
in Section [5] the deformations of various classical models. This will include the study of 
the Gaussian fluctuations away from the bulk for rather general Wigner and Wishart 
matrices, hence providing a new proof of the first part of [HI Theorem 1.1] and of [5l 
Theorem 3.1] but also a new generalisation to non-white ensembles. The study of the 
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1 .99 2.00 



(a) Case where 9 = 0.5 



(b) Case where — 1.5 



Figure 1. Comparison between the largest eigenvalues of a GUE 
matrix and those of the same matrix perturbed: the abscises of 
the vertical segments correspond to the largest eigenvalues of X, a GUE 
matrix with size 2.10 3 (under the dotted line) or to those of X n = X + 
diag(#, 0, . . . , 0) (above the dotted line). In the left picture, 9 = 0.5 < 9 = 1 
and as predicted, Ai ~ b = 2, whereas in the right one, 9 — 1.5 > 9, which 
indeed implies that Ai ~ po = 9 + \ = 2.17 and A2 ~ b. Moreover, in the 
left picture, we have, for all i, Aj ~ Aj, with some deviations 

W — W deviation of Aj from its limit 2. 

In the same way, in the right picture, i, X i+ i f» Aj, with some deviations 

|Aj + i — Aj| <C deviation of Aj from its limit 2. 

At last, here, in the right picture, we have Ai ~ 2.167, which gives 

Vn(\\-p e ) 



C0 



0.040, reasonable value for a standard Gaussian variable. 



eigenvalues that stick to the bulk requires a finer control on the eigenvalues of X n in the 
vicinity of the edges of the bulk, which we prove for random matrices such as Wigner and 
Wishart matrices with entries having a sub-exponential tail. This result complements [T8l 
Theorem 1.1], where the fluctuations of the largest eigenvalue of a non-Gaussian Wishart 
matrix perturbed by a delocalised but deterministic rank one perturbation was studied. 
One should remark that our result depends very little on the law v (only through its 
fourth moment in fact). 



Our approach is based upon a determinant computation (see Lemma 6.1 ), which shows 



that the eigenvalues of X n we are interested in are the solutions of the equation 



/„(*) := det ([Gl t (z)]l t=1 - diag^ 1 , • • ■ ,0) = 0, 



(3) 



with 



GUz) 



x„ 



"V>, (4) 
where (■, •) denotes the usual scalar product in C n . 

By the law of large numbers for i.i.d. vectors, by [101 Proposition 9.3] for uniformly 
distributed vectors or by applying Theorem 6.4 (with A n = (z — Xn)" 1 ), it is easy to see 
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that for any z outside the bulk, 

lim Gl t (z) = l s=t G, x (z) 

n— >oo 

and hence it is clear that one should expect the eigenvalues of X n outside of the bulk to 
converge to the solutions of G^ x (z) = 9J 1 if they exist. Studying the fluctuations of these 
eigenvalues amounts to analyse the behavior of the solutions of ^ around their limit. 
Such an approach was already developed in several papers (see e.g [7j or [I2])- However, 
to our knowledge, the model we consider, with a fixed deterministic matrix X n , was not 
yet studied and the fluctuations of the eigenvalues which stick to the bulk of X n was never 
achieved in such a generality. 

For the sake of clarity, throughout the paper, we will call "hypothesis" any hypothesis 
we need to make on the deterministic part of the model X n and "assumption" any hy- 
pothesis we need to make on the deformation R n . 

Moreover, because of concentration considerations that are developed in the Appendix 
of the paper, the proofs will be quite similar in the i.i.d. and orthonormalised models. 
Therefore, we will detail each proof in the i.i.d. model, which is simpler and then check 
that the argument is the same in the orthonormalised model or detail the slight changes 
to make in the proofs. 

Notations. For the sake of clarity, we recall here the main notations of the paper: 

• A" < < A™ are the eigenvalues of the deterministic matrix X n , 

• A" < < A™ are the eigenvalues of the perturbed matrix X n = X n + Y?i=i ^i v H u T , 

where r and the 0j's are independent of n and deterministic and the column vectors uf 
are random and defined above, 

• r e {0, . . . , r} is such that 9\ < ■ ■ ■ < 9 ro < < 9 ro+1 < ■ ■ ■ < 9 r , 

• for z out of the spectrum of X n , G™ t (z) = (u", [z — X n ) _1 ti"), 

• for z out of the support of fi, G^ x {z) = J -^dfj,x(x), 
m$ = 1— T > and 9 = T . < 0, 

lim^f, G Mx (z) — - lim zta Gfj. x (z) — ' 

• for any non null 9, 




G-l(l/9) if0e(-oo,0)U(0,+oc), 
a if 9 e [9,0), 

b if 9 e (0,9], 



• p + is the number of i's such that p g . > b, p_ is the number of z's such that p e . < a 
and cui < • • • < a q are the different values of the ^'s such that p e . ^ {a,b} (so that 
q < P- + p+, with equality in the particular case where the #j's are pairwise distinct), 

+ are the rescaled differences between the eigenvalues with limit out 
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of [a, b] and their limits: 

(y/HW-pet) if l<z<p_, 

7 H ~ 

l V^(K-( P - +P+)+ i - Pe r . {p _ +p+)+l ) if V- < * < V- + P+, 

• for any x such that x < a ot x > b, I x is the set of indices i such that = x, 

• for any j = 1, . . . , q, kj is the number of indices i such that Q{ = aj, i.e. kj = \I Pa . |. 

2. Almost sure convergence of the extreme eigenvalues 



For the sake of completeness, in this section, we prove Theorem L3 In fact, we shall 
even prove the more general following result. 



Theorem 2.1. Assume that Hypothesis LJ_ and Assumption \1.2\ are satisfied. 



Let us fix, independently of n, an integer i > 1 and V, a neighborhood of pg t if i < r 
and of a if i > r . Then A™ 6 V with overwhelming probability. 

The analogue result exists for largest eigenvalues: for any fixed integer i > and V, 
a neighborhood of pe r _ i if i < r — tq and of ' b if % > r — r^, A"_j G V with overwhelming 
probability. 



By Lemma |6.1| the eigenvalues of X n which are not in the spectrum of X n are the 
solutions of the equation 

det(M n (z)) = 0, 

with 

M n (z) = [Gi t (z)] r s>t=1 - dia g (0r\ . . . , e; 1 ), 

the functions G™ t (-) being defined in Q. For z out of the support of px, let us introduce 
the r x r matrix 

M{z) := diag(G MJf (z) - , G^z) - 9; 1 ). 



The key point, to prove Theorem 2.1 is the following lemma. For A = [Aij]lj =1 and rxr 
matrix, we set |y4|oo := sup^j \AiJ[. 



Lemma 2.2. Assume that Hypothesis \1.1\ and Assumption 1.2 are satisfied. For any 
5, e > 0, with overwhelming probability, 

sup \M{z)-M n {z)\ 00 <e. 

z, d(z,[a,b])>6 



In the case where the O^s are pairwise distinct, Theorem 2A_ follows directly from this 
lemma, because the z's such that det(M(z)) = are precisely the z's such that for some i, 
G flx (z) = J- and because close continuous functions on an interval have close zeros. The 
case where the ^'s are not pairwise distinct can then be deduced by an approximation 
procedure similar to the one of Section 6.2.3 of |lUj . 
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Proof of Lemma 2.2. The i.i.d. model. Fix R such that for all x £ [a — 5/2, b + 5/2] and 
z G C with \z\ > R, 



1 



z — X 



< 



Then since the support of px is contained in [a, b] and for n large enough, the eigenvalues 
of X n are all in [a — 5/2, b + 5/2], it suffices to prove that with overwhelming probability, 

sup \M{z)-M n {z)\ 00 <e. 

\z\<R,d(z,[a,b])>8 

Now, fix some z such that \z\ < R, d(z, [a, b]) > 5, and n large enough. By Proposition 
[z — X,^' 1 , whose operator norm is bounded by 25~ l , we find that for any 



6.2 with A 



e > 0, there exists c > such that 



P 



Gl t (z)- WTr((*-X 
n 



> 



5- 



l/2-e 



< 4e~ 



(5) 



It follows that there are c, rj > such that for all z such that \z\ < i?, <i(;z, [a, 6]) > 5, 

P(|Af(z)-Af n (^)| 00 >e/2)<e- cn \ 

As a consequence, since the number of z's such that \z\ < R and nz have integer real and 
imaginary parts has order n 2 , there is a constant C such that 

P(sup \M{z) -M n (z)\ OQ > e/2) < Cn 2 e~ cn \ 

z 

where the supremum is taken over complex numbers z = - + i-, with k, I £ Z, such that 
M < -R, ^l-^, [a, 6]) > 5. Now, note that for n large enough so that the eigenvalues of 
X n are all in [a — 5/2, b + 5/2], the Lipschitz norm for | • ^ on the set of ^'s such that 
d[z, [a,b]) > 5 of the function z \ — > M n (z) is less than 4. max s (=1 r || Therefore, 

again, with overwhelming probability z h 



by Proposition 6.2 again, with overwhelming probability z i — )■ M n (z) is -^-Lipschitz 
on this set. The function z i — > M(z) is 4-Lipschitz on this set, so, with overwhelming 
probability, 

sup \M n (z) - M(z)|oo < max \M n (z) - M(z)\ 00 + 85~ 2 n~ 1/2 , 

\z\<R,d(z,[a,b])>5 \z\<R,d(z,[a,b])>8 

nz&L+HL 

which insures that for n large enough, 

P f sup \M{z) -M n {z)\ oa > e\ < Cn 2 e~ cn \ 

\\z\<R,d(z,[a,b])>S J 

This concludes the proof for the i.i.d. model. 

The orthonormalised model can be treated similarly, by writing U n = W n G n with 
y/nW n a matrix converging almost surely to the identity by Proposition 6.3 □ 



3. Fluctuations of the eigenvalues away from the bulk 



3.1. Statement of the results. Let p + be the number of i's such that pg. > b and p_ 
be the number of i's such that pe i < a. In this section, we study the fluctuations of the 
eigenvalues of X n with limit out of the bulk, that is (A", . . . , A™_, A"_ p++1 , . . . , A^;). We 
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shall assume throughout this section that the spectral measure of X n converges to px 
faster than 1 / sfn. More precisely, 

Hypothesis 3.1. For all z £ {p ai , ■ ■ ■ ,Pa q }, \/n(G tln {z) — G^ x (zj) converges to 0. 

Our theorem deals with the limiting joint distribution of the variables 7™, . . . , 7p_ +J3+ , 
the rescaled differences between the eigenvalues with limit out of [a, b] and their limits: 

(V^(K-peJ if l<i<P- 

(V^(K-( P . +P+)+ i ~ P* r -( P _ +P+)+i ) lf P~ < 1 < V- +P+ 

Let us recall that for k > 1, GOE(/c) (resp. GUE(fe)) is the distribution of a k x 
symmetric (resp. Hermitian) random matrix [fl'tjli j=i such that the random variables 
{^M ! 1 < ^ < ^} U {^j ; 1 < i < j < k} (resp. {g^ ; 1 < % < k} U {v^^j) ! 1 < % < 
j < k} U {\/29f(<7i i3 -) ; 1 < i < j < fc}) are independent standard Gaussian variables. 

The limiting behaviour of the eigenvalues with limit outside the bulk will depend on 
the law v through the following quantity, called the fourth cumulant of v 

f x 4 dz/(x) — 3 in the real case, 
J |z| 4 dz/(^) — 2 in the complex case. 

Note that if v is Gaussian standard, then kAv) = 0. 



The definitions of the a/s and of the k/s have been given in Theorem 1.4 and recalled 
in the Notations gathered at the end of the introduction above. 



Theorem 3.2. Suppose that Assumption 1.2 holds with K±(y) = 0, as well as Hypotheses 
\1.1\ and 3J_. Then the law of 



converges to the law of (Xij, 1 < % < %)i<j<g, with Ajj i/ie z£/i largest eigenvalue of 
c a .Mj with (Mi, . . . , M q ) being independent matrices, Mj following the GUE(kj) (resp. 
GOE(kj)) distribution if v is supported on the complex plane (resp. the real line). The 
constant c a is given by 



T7 rlo-; — in the i.i.d. model, 

J {p a -x)- 2 d^ x {x) 



r d ^x( x ) 1 

^ (Pa-'-c)' 2 a 2 

k (/(P«-a:)- 2 d/ix(a;)) Z 



(6) 



in the orthonormalised model. 



When k±(v) ^ 0, we need a bit more than Hypothesis 3.1, namely 
Hypothesis 3.3. For all z £ M\[a, b], there is a finite number l(z) such that 
n S"=i(( z ~~ --^n) -1 )!? — ^ K z ) th e i-i-d. model, 



n—>oo 



Sr=i(((- 2 — X n ) l )a — - Tr((z — X n ) 1 )) 2 — > l(z) in the orthonormalised model. 
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We then have a similar result. 



Theorem 3.4. In the case when Assumption 1.2 holds with k±{v) ^ 0, under Hypotheses 
1.1, 3.1 and 3.3, Theorem 3.2 stays true, replacing the matrices c ai Mj by matrices c aj Mj + 



Dj where the Da 's are independent diagonal random matrices, independent of the Ma 's, 



and such that for all j, the diagonal entries of Dj 
variables, with variance —l{p aj )K^{v) / 'G' {p a .) . 



are independent centred real Gaussian 



3.2. Proof of Theorems |3.2| and |3.4[ We prove hereafter Theorem |3.2| and we will 

indicate briefly at the end of this section the minor changes to make to get Theorem 
3.4| The main ingredient will be a central limit theorem for quadratic forms, stated in 



Theorem 6.4 in the appendix. 



For i G {1, . . . q} and x G 1 
matrix with entries given by 



we denote by M n (i, x) the r x r (but no longer symmetric) 



[M n (i,x)] s y, 



We set p^x) := Pcti + ^. 

The first step of the proof will be to get the asymptotic behavior of M n (i,x). 



Lemma 3.5. Let i G {1, . . . q} and x G R be fixed. Under the hypotheses of Theorem 3.2 
M n (i,x) converges weakly, as n goes to infinity, to the matrix M.(i,x) with entries 

..... ,, ( G'„ x (fi a ,)(xt s=t -c n,. t ), ifs <=/ 



with (n s ,t)s,t=i,...,r & family of independent Gaussian variables with n SyS ~ A/"(0, 2) and 
n s ,t ~ A/"(0, 1) when s ^ t in the real case (resp. n SyS ~ jV(0, 1) and 9ft(n s ,t), ^s{n S)t ) ~ 
A/"(0, 1/2) and independent in the complex case). 



Proof. From ^ , we know that for s (fc I Pa . . 



Jim [M n (i,x)l t = [ G, x ( PcH ) - I) t s=t = fl - I ] l s=t . 



OsJ Vttj 9 S 



Let s G I Pa .. We write the decomposition 
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n [ «, {p l n {x) - X n )- l u n t ) - l s=4 -Tr((p^(x) - X n ) 

n 



-i- 



77 



"Tr(( Pc 



X r , 



The asymptotics of the first term is given by Theorem |6.4| with a variance given by 



lim -Tr{Mx)-X n 

n— >oo 77, 



G"u x {Pea) 



(9) 



As p ai is at distance of order one from the support of X n , we can expand xj^fn in 



M™t(i,x) to deduce that 



lim M^f(i, x) 



< x (p a ,)l s = t . 



Finally, by Hypothesis 3.1, we have 



lim M s f(i,x) = 0. 

n— >oo ' 



(10) 



Equations @, @, goj and (|TT} prove the lemma (using the fact that the distribution of 
the Gaussian variables n s>s and n Stt are symmetric). □ 

The next step is to study the behaviour of (M n (z, x)) xe K as a process on M. We will show 
in particular that the dependence in the parameter x is very simple. Let (n s t) s t=i...r De 



a family of Gaussian random variables as in Lemma |3.5| and define the random process 
•) from R to .M n (C) with [.M(i, x)],, t defined as in ^ (where we emphasize that 
(n s ,t)s,t=i,...,r do not depend on x). Then we have 

Lemma 3.6. Let i G {l,...q} be fixed. The random process (M n (i, x)) x£ ^ converges 
weakly, as n— >oo, to A4(i,-) in the sense of finite dimensional marginals. 



Proof. This is a direct application of Remark 6.5, as it is easy to check that for any 



x, x' e R, 



lim — Tr 

n— >oo 77, 



X 



-X, 



?? 



-i 



x V 
Pan H 7= — X T 

'n 



-r 



□ 



The last point to check is a result of asymptotic independence, from which the indepen- 
dence of the matrices M\,...,M q will be inherited. In fact, the matrices 
(M n (l, xi), . . . , M n (q, x q )) won't be asymptotically independent but their determinants 
will. 

Lemma 3.7. For any {x\, . . . ,x q ) G 1R 9 , the random variables 

det[M w (l,z 1 )],...,det[M n ( g> x,)] 
are asymptotically independent. 
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Proof. The key point is to show that, 

det[M n (i,x)]=det([M n (i,x)] s , teIpai ) J] (J- - j-) + (1) , (12) 

where the remaining term is uniformly small as x varies in any compact of 1R. 

Then, as the set of indices I Pai , ■ ■ ■ ,I Pa are disjoint, the submatrices involved in the 
main terms are independent in the i.i.d case and asymptotically independent in the or- 
thonormalised case. 



Let us now show (12). Firstly, note that by the convergence of M™ t {i } x) obtained in 
the proof of the Lemma 3.5, we have for all s, t £ {1, . . . , r} such that s ^ t or s £ I Pa , , 



for all k < 1 /2 

( G n .(n l (r)) - 

By the formula 



n K ( G™ t {p l n {x)) — l s =iTT- ) — > (convergence in probability). (13) 



/ n—too 



det [M»(i,x)] = n* £ s ^(cr) J] ( G^p^)) - 

aeS r s=l \ 

it suffices to prove that for any a £ S r such that for some Iq £ {1, . . . , r}\I Pa . , cr(i ) ^ i Q , 
* II ( G l°{s)(Pn( x )) - 1 *=*WfJ ^ (convergence in probability). (14) 

s=l ^ 



hi 
n 2 



It follows immediately from (13) since for any k < 1/2, in the above product, all the 
terms with index in I Pa . are of order at most n~ K , giving a contribution n~ ktK , and io is 
not in I Pa and satisfies cr(«o) ^ io, yielding another term of order at most n~ K . Hence, 
the other terms being bounded because p l n (x) stays bounded away from [a,b], the above 



product is at most of order n and so taking n £ ( gT^+ij i \) proves (14). □ 

Now as we have that, for i £ {1, . . . , q} and x £ 



det [M n (i, x)) = f n ( p ai + ^= ) 



jn 

we can deduce from the lemmata above the following 



Proposition 3.8. Under the hypothesis of Theorem 3.2, the random process 
converges weakly, as n goes to infinity to the random process 



G'^ip^detixl - c a M t ) II 



xeW l<i<q 
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in the sense of finite dimensional marginals, with the constants c ai and the joint distribu- 
tion of (Mi, . . . , M q ) as in the statement of Theorem 3.2, 



From there, the proof of Theorem 3.2 is straightforward. 
Proof. Let 

xi(i) < Vi(i) < x 2 (i) < y 2 (i) < ■ ■ ■ < y^i) 



be fixed. Since, by Theorem 21, for all e > 0, for n large enough, /„ vanishes exactly at 
P- + p + points in M\[a — e, b + s], we have that 



P 



xe(i) < l^-i k , f < ye(i), VI = I, . . . ,h, Vi = l,...q 

2-,m = l Krai-'- 



P 



fn Pa, + fn Pa t + 



<0, W =!,...,&*, Vz = l,...,g, 



-> P [det (y t (i)I - c ai Mi) det (x t (i)I - c a Mi) < 0, W = 1, . . . , k 4 , Vi = 1, . . . , q, } 

3 

= P [x e (i) < X it e < yi(i), V£ = 1, . . . , h, Vz = 1, . . . , q, ] 



□ 



To prove Theorem 3.4, the only substantial change to make is in the definition Q, in 
the case when s G I Pa . , we have to put 

[M(i,x)] s j := G' flx (p ai )(xl s=t - c Pa .n s>t ) - K 4 (v)l(p ai ). 

The convergence of [M n (i, x)] s t to [Ai(i, x)] st is again obtained by applying Theorem 



6.4 



4. The sticking eigenvalues 

4.1. Statement of the results. To study the fluctuations of the eigenvalues which stick 
to the bulk, we need a more precise information on the eigenvalues of X n in the vicinity of 
their extremes. More explicitly, we shall need the following additional hypothesis, which 
depends on a positive integer p and a real number a G (0, 1). Note that this hypothesis 
has two versions: Hypothesis 4.1[ p, a, a] is adapted to the study of the smallest eigenvalues 
(it is the version detailed below) and Hypothesis 4.1[ p, a, b] is adapted to the study of the 
largest eigenvalues (this version is only outlined below). 

Hypothesis 4.1. [p, a, a] There exists a sequence m n of positive integers tending to in- 
finity such that m n = 0(n a ), 



1 

lim inf — > 

n— >oo n 



1 1 

> -. 



i=m n + l P 

and there exist r] 2 > and 774 > ; so that for n large enough 



E 

i=m n +l 



(A? - AD 



< n 



2-m 



and 



E 



i=m n + 



1 (K ~ K) 4 



< n 4 - m 



(15) 

(16) 
(17) 
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Hypothesis 4.1 [p, a, b) is the same hypothesis where we replace A™ — A" by A™_ p+1 — 



A"_j +1; and (1151) becomes 



1 - 
lim sup — 

n— ¥oo W> . 



i+1 A n-p+l 



A n-i+l 



< 



For many matrix models, the behaviors of largest and smallest eigenvalues are similar, 



and Hypothesis 4.1 [p, a, a] is satisfied if and only if Hypothesis 4.1| [p, a, b] is satisfied. In 
such cases, we shall simply say that Hypothesis 4.1 [p, a] is satisfied. 



For rank one perturbations and in the i.i.d. model, we will only require the two first 



conditions (15) and (16) whereas for higher rank perturbations, we will need in addition 



(17) to control the off-diagonal terms of the determinant. 



Moreover, we shall not study the critical case where for some i, 9j G {9, 9}. 
Assumption 4.2. For all i, 9i ^ 9 and 9i ^ 9. 



In fact, Assumption 



4.2 



can be weakened into: for all i, 9^ ^ 9 (resp. 9i ^ 9) if we only 
study the smallest (resp. largest) eigenvalues. 

The fact that the eigenvalues of the non-perturbed matrix are sufficiently spread at the 
edges to insure the above hypothesis allow the eigenvalues of the perturbed matrix to be 
very close to them, as stated in the following theorem. 

Theorem 4.3. Let I a = {i G [l,r] : pg i = a} = [p_ + l,r ] (resp. 4 = {i 6 [l,r] : 
Pq. = b} = [r + 1, r — p + ]) be the set of indices corresponding to the eigenvalues A™ (resp. 
A"_ r+iy ) converging to the lower (resp. upper) bound of the support of px- Let us suppose 
Hypothesis 1.1, Hypothesis 1^.1 [r, a, a] (resp. Hypothesis 1^.1 [r, a, b]) and Assumptions 



1.2 and 4.2 



to hold. Then for any a' > a, we have, for all i G I a (resp. i G lb), 



min |A™ - A£| < n 

l<k<i+r—ro 



-l+a' 



(resp. 

with overwhelming probability. 



min |A"_ r+i - A£| < n 

n—r+i—ro<k<n 



-l+a' 



Moreover, in the case where the perturbation has rank one, we can locate exactly in 
the neighborhood of which eigenvalues of the non-perturbed matrix the eigenvalues of the 
perturbed matrix lie. 

We state hereafter the result for the smallest eigenvalues, but of course a similar state- 
ment holds for the largest ones. 

Theorem 4.4. Let (A")j>i be the eigenvalues of X n + 9u\u\, with 9 < 0. Then, under 
Assumption 1.2 and Hypothesis \1.1\ if (15) and (16) in Hypothesis 4-1 \p, ol, a] hold for 
some a G (0, 1) and a positive integer p, then for any a' > a, we have 

(i) if 9 < 9, A" converges to pg < a whereas n 1_Q '(A™ +1 — A")i<j< p _i vanishes in prob- 
ability as n goes to infinity, 
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(ii) if 8 £ (0,0), n 1_Q/ (A™ — A")i<j< p vanishes in probability as n goes to infinity, 



(iii) if, instead of (15) and (16) in Hypothesis 4-1 [p,ce,a], one supposes (15) and 
(16) in Hypothesis 4-1 \p,ct,b] to ZioZd, i/ien n 1_Q '(A™„ i — A"_Jo<i<p vanishes in 



probability as n goes to infinity. 

Theorem 4.5. Consider the i.i.d. model and let (Af)t>i fre i/ie eigenvalues of X n + 
X^=i 6i u i u i ■ Let p_ (resp. p + ) be the number of indices i so that pg. < a (resp. pg i > 
b). We assume that Assumptions 1.2 and 4-2, Hypothesis \1.1\ and (15) and (16) in 
Hypotheses 4-1 [p, a] and [q, a, b] hold for some a £ (0, 1) and integers p, q. Then, for 
all a' > a, for all fixed 1 < i < p — (p_ + r) and < j < p — (p + + r), 



n 1 ~ a (Xp_ +i — A") and ^(A^^ 
&otA vanish in probability as n goes to infinity. 



Ki-i) 



Note that if p — (p_ + r) < (resp. if p — (p + + r) < 0), then the statement of the 
theorem is empty as far as z's (resp. j's) are concerned. The same convention is made 
throughout the proof. 



4.2. Proofs. Let us first prove Theorem |4.3| Let us choose Iq £ I a and study the be- 
haviour of A" o (the case of the largest eigenvalues can be treated similarly). We assume 
throughout the section that Hypotheses 1 1 . 1 1 |4~T| [r, a, a] and Assumptions 1.2 and 4.2 
satisfied. We also fix a' > a. 



are 



We know, by Lemma 6.1 that the eigenvalues of X n which are not eigenvalues of X n 
are the z's such that 

det(M n (z)) = 0, (18) 

where 

M n {z) = [Gl t (z)]l t=1 -diag^r 1 ,-..,^ 1 ) (19) 

and for all s,t, 

Gl t (z) = (u:,(z-X n y 1 u?). 
Recall that by Weyl's interlacing inequalities (see [H Th. A. 7]) 

XI < K 



x i — ^i +r-r - 



Let C be a fixed constant such that maxi<j< p _ 
that 



< £ < a. By Theorem 2.1, we know 



Lemma 4.6. With overwhelming probability, A™ > £ 



We want to show that (|18|) is not possible on 



K'-=\ze [C, A™ +r _ r J ; mm 

I l<k<io+r—ro 



\z - At I > n 



-l+o' 



The following lemma deals with the asymptotic behaviour of the off-diagonal terms of 



the matrix M n (z) of (19). 
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Lemma 4.7. For s ^ t and k > small enough, 

sup \Gl t (z)\ <n~ K 

z£fi n 



(20) 



with overwhelming probability. 

The following lemma deals with the asymptotic behaviour of the diagonal terms of the 
matrix M n (z) of Q. 

Lemma 4.8. For all s = 1, . . . , r, for all 5 > 0, any 5 > 0, 

1 



sup 



< S 



(21) 



with overwhelming probability. 

Let us assume these lemmas proven for a while and complete the proof of Theorem 



4.3| By these two lemmas, for z 6 fi n , we find by expanding the determinant that with 
overwhelming probability, 



det(M n (z)) = H(Gl 



8=1 



+ 0(n~ K ) 



(22) 



where the 0(n K ) is uniform on z G Q n . Indeed, in the second term of the right hand 
side of 



8=1 ^ l ' a£S r \{Id} S=l 

each diagonal term is bounded and each non diagonal term is 0(n~ K ). 



i S =a(s)9 s X ), 



Since for all i, 9i ^ 9, (22) and Lemma 4.8 allow to assert that with overwhelming 
probability, for all z G O n , det(M n (z)) ^ 0. It completes the proof of the theorem. □ 

We finally prove the two last lemmas. 



Proof of Lemma ^.7. Let us consider z G fi n (z might depend on n, but for nota- 
tional brevity, we omit to denote it by z n ). We treat simultaneously the orthonormalised 
model and the i.i.d. model (in the i.i.d. model, one just takes W n = I and replaces 
||(G ?n (Vy n ) T ) s ||2 by y/n in the proof below). Observe that if we write X n = 0*D n O with 
D n = (A", . . . , A") and O a unitary or orthogonal matrix, 

guz) = (u:,(z-x n yx) 



£ (Oti?),(Oti?), 



1=1 



z-\? 



The first step is to show that for any e > 0, with overwhelming probability, 



max | {Ou r & 

l<n,s<r 



< n 2" 



(23) 



Indeed, with 0\ the Zth row vector of O and using the notations of Section 6.2 



\\(G n (W 



n\T\ 



)s\\2 



t=l 
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But g i — y (Oi,g™) is Lipschitz for the Euclidean norm with constant one. Hence, by 
concentration inequality due to the log-Sobolev hypothesis (see e.g. [U section 4.4]), 
there exists c > such that for all 5 > 0, 



F(\(O h g:)\>6)<4e 



cS 2 



so that 



P ( max \{O h g n s ) \ > n e ) < 4n 4 e 

Kn,s<r 



4„— cn 2 



From Proposition 6.3, we know that with overwhelming probability, ||(G n (W /n ) ) s ||2 is 



bounded below by ^Jnn 6 and the entries of W n are of order one. This gives therefore 



(23). 



We now make the following decomposition 



GUz) 



E 

i=i 
\ 



z - X 



l=m n +l 



z-Xf 



■■=A n (z) 



— V 

=B„(z) 



Note that as |(Ou™);|, 1 < Z < m n , are smaller than n 2+ e by (23), for any e' > 0, with 
overwhelming probability, we have, uniformly on z £ Q n , 

\A n (z)\ < m n n 1_a V 1+2e ' = 0(n a - a>+2e ') 

We choose < e' < (a' — «)/4 and now study B n (z) which can be written 

B n (z) = (u:,P(z-X n )- 1 Pul l ) 

with P the orthogonal projection onto the linear span of the eigenvectors of X n corre- 



sponding to the eigenvalues A^ +1 , . . . , A™. By the second point in Proposition 6.2, with 
z £ Q n , for all s ^ t, 

< 4e _c<52 + 4e~ cmin ( K ' K2 \ 



Moreover, by Hypothesis 4.1 for n large enough, for all z £ £l n , 

Ti{P{z - X n )~ 2 ) < n 2 ~ m and Tr(P(z - X n )" 4 ) < n 4 "" 4 . 
We deduce that there is C, r\ > such that for all z £ fi n , 

1 



P 



n 



■i£,P{z-X n )-*P$) 



> n s < Ce 



A similar control is verified for s = t since we have, by Proposition |6.2 
1 (g s , P(z - X n )- x Pg.) - -Tr (P(z - X^ 1 ) 



P 



n 



71 



> S) <4e 



cmin{5 2 n''2 ,5^2/2} 



(24) 



(25) 



whereas Hypothesis 4.1 insures that the term ^-Tr(P(z — X n ) _1 ) is bounded uniformly on 
Q n . Thus, up to a change of the constants C and rj, there is a constant M such that for 
all z £ Q n , 

' P(z-X n )- l Pg s ) 



P 



n 



> M \ < Ce~ n \ 
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Therefore, with Proposition 6.3 and developing the vectors w"'s as the normalised column 
vectors of G n (W n ) T , we conclude that, up to a change of the constants C and rj, for all 

z G ^715 

P (\B n (z)\ > rT 22 ^) < Ce~ n \ (26) 

Hence, we have proved that there exists k > 0, C and rj > so that for all z G Q n , 

F(\Gl t (z)\ > n~ K ) < Ce~ n \ 

We finally obtain this control uniformly on z G Q n by noticing that z— >G™ t (z) is Lipschitz 
on Q n , with constant bounded by (min|z — Aj|)~ 2 < n 2 ~ 2a ' . Thus, if we take a grid 
( z k )o<fc<cn 2 of Q n with mesh < n~ 2+2a ~ K (there are about n 2 such zj^s) we have 

sup |G? (z)| < max jG? it (*j?)| + n"\ 
Since there are at most cn 2 such k and n 2 possible z, j, we conclude that 



P (sup \Gl t (z)\ > 2n- K ) < c 2 n 4 Ce- nV 



which completes the proof. □ 

Proof of Lemma \4-8[ We shall use the decomposition 

GUz) = «, P(z - X n y l Pv%) + «, (1 -P)(z- X n y\l - P)<), (27) 

with P as above the orthogonal projection onto the linear span of the eigenvectors of X n 
corresponding to the eigenvalues A^ n+1 , . . . , A™, and then prove that for z G Q n , 



whereas 



(u^p{z-x n )- l Pu n s )* 1 -, 



i 



(«?,(! -PXz-An^Cl-PK) < — — 11(1 - 



n\\2 
s 1 1 2 



mini< fc < mn \z-X\ 

v v ' Sin" 1 rank(l-P) 

<n 1 - a ' 



« n- a 'm n = 0(n a - a ') = o(l). 

Let us now give a formal proof. Again, we first prove the estimate for a fixed z G O n , 
the uniform estimate on z being obtained by a grid argument as in the previous proof (a 
key point being that the constants C and r\ of the definition of overwhelming 'probability 
are independent of the choice of z G fl n ). 



First, observe that (|15|) implies that for any sequence e n tending to zero 

1 _l 1 

n 2-^ y — \n ft 



lim sup 

n^oo a-e„<z<\V, 



n ^ z-Xf 9 

i=m n + l 



0. (2* 



Indeed, for all e > 0, for n such that A^ and a — e n are both > a — e, we have, for all 
z G [a - e n ,Xp], 

l A ii" ii" i 

n 2-" A" - A? ~ n ^ ^ - A" ~ n ^ a - e - A™ ' 

i=m n +l f 1 i=m n + l 1 i=m n +l 1 
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so that ( 15 ) and 



1 n 

lim lim sup — > — 

40 n^oo n . ^ a 

i=m n +l 



X? 



1 

6 



imply (|28j). 

So let us consider z G Q n (z might depend on n, but for notational brevity, we omit to 
denote it by z n ). By the inequality \z — A£| > n~ 1+Q for all 1 < k < m n and (27), we 
have 



\GUz) - (u:,P(z - X n )- l Pu:)\ < n l ~« 11(1 - P)u n s \\l (29) 



But as in the previous proof, we have 



||(G"(W 



n\T\ 



ls\\2 



z — / i n 



Lv=l 



with, by (24), the off diagonal terms t ^ v of order n~ mAm ' 8 with overwhelming prob- 
ability, whereas the diagonal terms are close to -Tr(P(2 — X n ) _1 ) with overwhelming 
probability by (25). Hence, we deduce with Proposition 6.2 that for any 5 > 0, 



(u:,p(z-x n )- i p^) 



1 



77. 



Tr(P((z - X n )~')) 



< 5 



with overwhelming probability. Hence, by (28), for any 5 > 

1 



[u n sl P{z-X n )^Pu n s ) 

with overwhelming probability. On the other hand 

1 







< 5 



(30) 



P)u 



n || 2 
2 



ll(^(w-F) s ||^ ( 



6.3[ the denominator is of order n with overwhelming probability, whereas 



By Proposition 
by Proposition 

overwhelming probability. As W n is bounded by Proposition |6.3| we conclude that 



6.2, the numerator is of order m n + rfJm n (since Tr(l — P) = m n ) with 



PK\\i<2- 



n 



(31) 



with overwhelming probability. Putting together Equations (29), (30) and (31), we have 
proved that for any z G Q n , any 5 > 0, 

1 

6 



< 5 



with overwhelming probability, the constants C and r\ of the definition of overwhelming 
'probability being independent of the choice of z G VL n We do not detail the grid argument 
used to get a control uniform on z because this argument is similar to what we did in the 
proof of the previous lemma. □ 



Proof of Theorem 4^ in the one dimensional case, the eigenvalues of X n which do not 
belong to the spectrum of X n are the zeroes of 



f n {z) = -(g,{z-X n ) 1 g)-e n (g)- 
n v 



(32) 



EXTREME EIGENVALUES OF DEFORMED RANDOM MATRICES 



21 



with e n (g) = 1 or ||<7|||/ti according to the model we are considering. A straightforward 
study of the function f n tells us that the eigenvalues of X n are distinct from those of X n 
as soon as X n has no multiple eigenvalue and 

(matrix of the eigenvectors of X n )* x g 

has no null entry, which we can always assume up to modify X n and g so slightly that 
the fluctuations of the eigenvalues are not affected. We do not detail these arguments but 
the reader can refer to Lemmas 9.3, 9.4 and 11.2 of [H] for a full proof in the finite rank 
case. 



Therefore, (|32|) characterises all the eigenvalues of X n . Moreover, by Weyl's interlacing 

\ n l <\ n 1 <\ n 2 <\ n 2 <---<\ n n <K- 



properties, for 9 < 0, 



Theorems 2.1 and 4.3 thus already settle the study of A" which either goes to pg or is 
at distance 0(n~ 1+a ) of A™ depending on the strength of 9. We consider a' > a and 
i G {2, . . . ,p} and define 



A n := 



+ n- 1+a \ A? 



n 



-l+a' 



Note first that if A n is empty, then the eigenvalue of X n which lies between Af_ x and A™ is 
within n~ l+a ' to both A™_ x and A™, so we have nothing to prove. Now, we want to prove 
that f n does not vanish on A n and that according to the sign of \ — 4, it vanishes on one 
side or the other of A n in ]A£_ 1; A™[. This will prove (i) and (ii) of the theorem. Part (iii) 
can be proved in the same way, proving that with overwhelming probability, /„ does not 
vanish in ] X^-i-i + n~ 1+a ' , A™_; - n- 1+a ' [. 



The proof of this fact will follow the same lines as the proof of Lemma 4.8 and we 
recall that P was defined above as the orthogonal projection onto the linear span of the 
eigenvectors of X n corresponding to the eigenvalues A^ +1 , . . . , A™. Then, exactly as for 



pO) , we can show that for all 5 > 0, 

1 , 



sup 



n 



P(z-X n )- L Pg) 



< S 



with overwhelming probability. Moreover, for any z G A n , for any j 



, m n , we have 



- X]\ > mm{z - A"_ l5 A™ - z] > n 



-l+a' 



so that 



sup 

z£A n 



n 



(g,(l-P)(z-X n )- 1 (l-P)g) 



<n- a '(g,(l-P)g). 



By Proposition 6.2[ we deduce that for any e > 

1 



sup 



n 



,{\-P){ z -X n )- 1 {l-P)g) 



< n e n 



m r , 



with overwhelming probability. We choose e in such a way that the latter right hand side 
goes to zero. Therefore, we know that uniformly on A n , 

1 1 



fn{z) 



9 9 



+ o(l) 
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with overwhelming probability. Since for all n, f n is decreasing, going to +00 (resp. 
—00) as z goes to any A™_ x on the right (resp. A™ on the left), it follows that according 
to the sign of 4 — 4, the zero of f n in ]A"_ 1; A™[ is either in ]A™_ 1 ,A™_ 1 + n~ l+a '[ or in 
]A™ - n- 1+a> , A™[. □ 
Proof of Theorem \4-5[ For each £ = 0, . . . , r, let us define 



i=i 

and denote its eigenvalues by < < An ■ We also define 

p_ := ${i = I,. . . p 6l < a}, 
pf ■= t{i = l,...,l;po i >b}. 

(£) (£) ' — CO 

p_ and p + are respectively the numbers of eigenvalues of X n with limit < a and > b. 
We also set 

Of course, as before, the eigenvalues of Xn ^ are the zeros of fn\ 
Let us also choose ( a < a and Cb> b such that 

C a > max{p 9i ; p 6i < a} and ( b < mm{p dz ; p e . > b}. 

First, as in the proof of Theorem |4.4[ up to small perturbations, one can suppose that 
for all £ = 0, . . . , r, the eigenvalues of X n are pairwise distinct and for all £ = 1, . . . , r, 

~{£) —(£-1) 

the eigenvalues of X n are distinct from those of X n 
Now, let us state a few facts: 

— (0 

(a) For all £, there is a constant M such that the extreme eigenvalues of X n are in 
[—M,M] with overwhelming probability (this follows from Theorem 2.1). 

(b) Moreover, for each £, for each i < k, by Weyl's interlacing inequalities, 

o- 1 • 1 



A i+i _ \-i \ ~ A k 



which implies, by induction over £, that X n satisfies the first part Hypothesis and 



— (£) ~(£-l) 

We only consider the i.i.d. model, so each X n can be deduced from X n by adding 



( 15 ) and ( 16 ) in Hypotheses |4.1| \p — £, a, a] and [q — £, a, b] 

We only consider the i.i.d. model, so 
an independent rank one perturbation. 

~ CO 

In the case where all the ^'s are in [9,9], also the extreme eigenvalues of X n stick 
to the bulk and therefore the full hypothesis LI holds at each step. In this case we can 

y to prove the theorem. The appearance of spikes is 
insures that for all I, the eigenvalues of Xn are close 



simply apply Theorem L4 inductive 
in fact not a problem as Theorem 4.3 



to the eigenvalues of X n simultaneously with overwhelming probability, whereas Weyl's 
interlacing properties and as in the previous proof discussions on the sign of the functions 
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/„ allows to localise in the neighborhood of which eigenvalues of X n the eigenvalues of 
A-;," lie. 

Let us d etai l a bit this argument. By using (a) and (b) above and following the proof 
of Lemma 4.8l one can easily prove that for all I = 1, . . . , r, for any 
. . . , q — £), for any 5 > 0, for 

,-l+a'r 



4.8 

(resp. j = p+ 



P- , 



P 



n 



(resp. n n :=]A^ 3 . + min{A W 
with overwhelming probability, 



sup 



(z-X, 



(£-1) 



1,0,} 



n 



-l+a'l 



(resp. sup 



< 5 



<*)■ 



Let us now fix £ G {l,...,r} and compare the eigenvalues of X n to the ones of X n 
We suppose for example that 9e > 0. 
Then by Weyl's inequalities, we have 

Xt 1] < A? < A 



(33) 
(34) 

(35) 

(36) 
(i-i) 



<xtP<x^l<x^<\^. 



v n-l 



Let us first consider the smallest eigenvalues. Under the overwhelming event (35), 
fn < on any interval Q n as defined in (33). So, since /„ is decreasing and 
vanishes exactly once on ]Ajf~[ , A^ ^f, its zero Xf\ is within n~ 1+a ' from A^ 1 ^. 



Let us now consider the largest eigenvalues. Under the overwhelming event (|36|) 
ame sig 
from A 



fn has the same sign as = — j- on any interval f2 n as defined in (34), so A„i~is 



within n 



, ii 9i < 9 and A„_^- is within n 



from A^T/la if 9 X > 9. 



To conclude, up to n~ 1+a ' errors, each perturbation by a positive rank one matrix 
9eii™u™ does move the smallest eigenvalues and translates each largest one to the following 
eigenvalue if 9\ > 9 and does not move the largest eigenvalues if 9i < 9. Of course, the 
analogue result holds for perturbations by negative rank one matrices. 

The theorem follows. □ 



5. Application to classical models of matrices 



Our goal in this section is to show that if X n belongs to some classical ensembles of 
matrices, the extreme eigenvalues of perturbations of such matrices have their asymptotics 
obeying to Theorems 2.1 3.2 and 4J3 For that, a crucial step will be the following 
statement. If (X n ) is a sequence of random matrices, we say that it satisfies an hypothesis 
H in probability if the probability that X n satisfies H converges to one as n goes to infinity 
(for example, if H states a convergence to a limit £, "if in probability" is the convergence 
in probability to £). 
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Theorem 5.1. Let (X n ) be a sequence of random matrices independent of the u™ 's. Under 



Assumption 1.2 



(1) If Hypothesis 

(2) If k±(v) = and Hypotheses 



1.1 holds in probability, Theorem 2.1 holds 



1.1 and 3.1 hold in probability, Theorem 3.2 holds 



1.1 and 3.3 hold in probability, Theorem 3.4 holds 



7/k 4 (z/) and Hypotheses 
(3) Under Assumption 4-2, if Hypotheses 1.1 and 4-1 hold in probability, Theorem 4-3 
holds "with probability converging to one" instead of "with overwhelming probabil- 
ity"; Theorems \4-4\ and Corollary hold. 



This result follows from the results with deterministic sequences of matrices X n . In- 
deed, to prove that a sequence converges to a limit £ in a metric space, it suffices to prove 
that any of its subsequences has a subsequence converging to I. If the convergences of the 
hypotheses hold in probability, then from any subsequence, one can extract a subsequence 
for which they hold almost surely. Then up to a conditioning by the u-algebra generated 
by the XJs, the hypotheses of the various theorems hold. 



The remaining of this section is devoted to showing that such results hold if X n , inde- 
pendent of (w")i<j< r , is a Wigner or a Wishart matrix or a random matrix which law has 
density proportional to e~ TrV for a certain potential V. In each case, we have to check 
that the hypotheses hold in probability. 



5.1. Wigner matrices. Let /ii be a centred distribution on R (respectively on C) and 
\ii be a centred distribution on IR, both having a finite fourth moment (in the case where 
fii is not supported on the real line, we assume that the real and imaginary part are 
independent). We define a 2 = J z£C |z| 2 d/zi(,2). 

Let {xi,j)i,j>i be an infinite Hermitian random matrix which entries are independent 
up to the condition such that the distributed according to /z 2 and the 

Xij's {i 7^ j) are distributed according to H\. We take X n = 4= [xij]^- =l , which is said 
to be a Wigner matrix. For certain results, we will also need an additional hypothesis, 
which we present here: 

Hypothesis 5.2. The probability measures fii and /i2 have a sub- exponential decay, that 
is there exists positive constants C,C such that if X is distributed according to /ii or fi 2 , 
for all t > C, 

H\ x \ > fC ) < e ~ l - 

Moreover, fi\ and \ii are symmetric. 



The following Proposition generalises some results of [36l HI CGI [13] which study the 
effect of a finite rank perturbation on a non-Gaussian Wigner matrix. In particular, it 
includes the study of the eigenvalues which stick to the bulk. 



The 



Proposition 5.3. Let X n be a Wigner matrix. Assume that Assu mption 1.2 holds, 
limits of the extreme eigenvalues of X n are given by Theorem 2. 1 and the fluctuations of 
the ones which limits are out of [— 2cr, 2<r] are given by Theorem 3.2, where the parameters 
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a, b, p e , c a are given by the following formulas : b = —a = 2a, 

if\6\ > a, 
p e := { 2a ifO < 9 < a, 



and 




if-a<9 < ; 

\/a 2 — a 2 in the i.i.d. model, 

TVa ~ a in the orthonormalized model. 



Assume moreover that, for all i, 9{ ^ {—a, a} and Hypothesis 5.2 holds. If the pertur- 
bation has rank one, we have the following precise description of the fluctuations of the 
sticking eigenvalues : 

• If9>a (resp. 9 < -a), for all p > 2, n 2/3 (A"_ p+1 - 2a) (resp. n 2 / 3 (A™ + 2a)) 
converges in law to the p — 1th Tracy Widom law. 

• If < 9 < a (resp. -a < 9 < 0), for all p > 1, n 2/3 (\™_ p+1 - 2a) (resp. 
n 2//3 (A™ + 2a) ) converges in law to the pth Tracy Widom law. 

If the perturbation is rank more than one and Assumption \4-2\ holds, the extreme eigen- 
values of X n are at distance less than n~ 1+e for any e > to the extreme eigenvalues of 
X n , which have Tracy- Widom fluctuations. We can localize exactly near which eigenvalue 
of X n they lie by using Theorem \4-5\ in the i.i.d model. 

Remark 5.4. All the Tracy-Widom laws involved in the statement of the proposition 
above, are the ones corresponding respectively to the GOE if p,\ is supported on IR and to 
the GUE if pi is supported on C. 



According to Theorem |5.1[ it suffices to verify that the hypotheses hold in probability 
for (X n ) n >i. We study separately the eigenvalues which stick to the bulk and those which 
deviate from the bulk. 

•Deviating eigenvalues. 

If X n is a Wigner matrix (that is, with our terminology, with entries having a finite 
fourth moment), the fact that X n satisfies Hypothesis in probability is a well known 
result (see for example jH Th. 5.2]) for p,x the semicircle law with support [—2a, 2a}. 
The formulas for pg and c a can be checked with the well known formula |pj[ Sect. 2.4]: 

Wz e R\[-2a,2a], G, x {z) = z '^(zWz 2 -^ ^ 

Moreover, [5j Th. 1.1] shows that Tr(f(X n )) - nj f(x)da( x) converges in law to a 
Gaussian distribution for any function / which is analytic in a neighborhood of [—2a, 2a]. 
For any fixed z ^ [—2a, 2a], applied for f(t) = -^r, we get that niG^z) — G ilx {z)) 
converges in law to a Gaussian distribution, hence y/n(G^ n (z) — G^ x {z)) converges in 



probability to zero, so that Hypothesis 34 holds in probability. 
• Sticking Eigenvalues. 
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We now assume moreover that the laws of the entries satisfy Hypothesis 5.2 In order 
to lighten the notation, we shall now suppose that a — 1. Let us first recall that by 
[HI [39], the extreme eigenvalues of the non-perturbed matrix X n , once re-centred and 
renormalised by n 2 ^ 3 , converge to the Tracy- Widom law (which depends on whether the 
entries are complex or real). We need to verify that Hypothesis 4.1 p,a] for any finite p 
and an a < 1/3 is fulfilled in probability. By jH], the spacing between the two smallest 
eigenvalues of X n is of order greater than n -7 for 7 > 2/3 with probability going to one 
and therefore, by the inequality 



E 



t=m n + 



n\l—k 



E 

i=m n +l 



X? - XI 



(k = 2 or 4), 



it is sufficient to prove the first point of Hypothesis 4.1 p,a]. We shall prove it by replacing 
first the smallest eigenvalue by the edge —2 thanks to a lemma that Benjamin Schlein 
[4*0] kindly communicated to us. We will then prove that the sum of the inverse of the 
distance of the eigenvalues to the edge indeed converges to the announced limit, thanks 
to both Soshnikov paper |4T] (for sub-Gaussian tails) or [39] (for finite moments), and 
Tao and Vu article 



Lemma 5.5 (B. Schlein). Suppose the entries of X n have a uniform sub- exponential tail. 
Then for all 5 > 0, for all integer number p, 



lim P 

n— >oo 



I " 1 1 " _l 



n A™ 
j= P +i 3 



n A? + 2 
j=p+i J 



> S 



0. 



Proof. We write 

n 

- y — - 



n X r - 
j= P +i 3 



A I 



n 1 

-Y — 

n ^ XI- 



j=p+i 



A™ + 2 



A™ + 2 " 



n 



3=P~- 



1 (A?->?)(Ay + 



Hence for any K\ > 
P 



[1 A 1 1 ^ 1 

\ j=p+l 3 P j= p+ l 3 



> s 



< 



\x; + 2\ > Km 




-2/3N 



V i 



> 5 and |A™ + 2| < K x n 



Now, for any K 2 > K±, on the event {|Ap + 2| < Kin 2//3 }, for any k > 0, we have 



-2/3 



(3? 



Ky 

n 5/3 



E 



p 

+00 



£ x I(a?-W + 2)I 



< 



Ki ¥Z Af n [2K 2 n- 2 / 3 + £n- K , 2K 2 n~ 2 / 3 + (£ + l)n~ K } 



£=0 



{K 2 n~ 2 l 3 + £n- K f 



5/3 



n 



E 



1 



j=p+l lv 3 



X j +2<2K 2 n- 2 / 3 

|(A«-A«)(A« + 2)| 



(39) 



where A4K &] := tl{* i — 2 + a < A™ < —2 + 6}. Note that, from the upper bound on the 
density of eigenvalues in microscopic intervals, due to [151 Theorem 4.6], we know that 
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for any k < 1, there is a constant M independent of n so that for all £ > 1 

E(A^„[2ir 2 n- 2/3 + 2K 2 n~ 2/3 + (£ + 1)t2~ k ]) < Mn 1 ~ K . (40) 



Let us fix k £ (g, 1). It follows that the first term of the r.h.s. of (39) can be estimated 
by 



P 



Ki ^Af n [2K 2 n- 2 / 3 + £n- K ,2K 2 n- 2 / 3 + (£+l)n- K } 5 



n 



5/3 



E 



(K 2 n~ 2 / 3 + In-*) 2 



> 



< g E(Ar n [2K 2 n - 2 / 3 + £n~ K , 2K 2 n~ 2 l 3 + [£ + l)n- K ]) 



2MKi 1 



(K 2 n~ 2 / 3 + £n~ K f 



Sn 2/3 u k Z_, (K 2 n~ 2 / 3 + £n—y 



< 



2MK l 



2MK l 



6n 2 /z n*(K 2 n- 2 / 3 ) 2 ' 8n l J (t + K 2 n~l) 



+oo 



dt 



2MK X 2MK X 
~ 5K 2 n^ 2 / 3 + 5K 2 ' 



(41) 



Let us now estimate the second term of the r.h.s. of (39). For any positive integer K%, 
we have 

K\ l|A^+2|<2if 2 n- 2 /3 5 



< P (K(-oo, 2K 2 n- 2 /' 3 ] > K 3 ) + P ( 



1 



< P(A£ < -2 + 2K 2 n- 2/3 ) + P min |A™ + 2|< 



6 

n^mm^^K, |(A? - A") (Ay + 2)| " 2 



v 7 ^ 



+p(Ia;-a; +1 |< 



71 



From (|38]), g9j, gy and g2j, we conclude that 

1*1 1 A 1 

P - > > > 5 

n ^ A? - A? n ^ A™ + 2 



(42) 



< PflAJ + 2| > K^ 2 / 3 ) + + P (\ Ka < -2 + 2iT 2 n- 2 / 3 ) 

+P( min |A™ + 2|<- ^ ] + — + P I |A 2 - A"| < - A 



\^<j<K 3 



V6 



SK 2 



for arbitrary < K\ < K s and K 3 > 1. Taking the limit n — > 00, the last two terms 
disappear, because by [121 Th. 1.16], the distribution of the smallest K 3 eigenvalues lives 
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on scales of order n 2 1 3 ^> n 5 / 6 . Therefore, 



lim P 

n— >oo 



lA 1 lA 1 



i=2 



n <7-Z A" + 2 

j=2 J 



> 5 



< lim P(|A™ + 2| > K in - 2/3 ) 



2MK 1 
SKo 



+ lim P (X Ka < -2 + 2K 2 n- 2/3 ) 



still for any < K\ < K% and Ks > 1. Now, note that for K\ large enough, the first term 
can be made as small as we want. Then, keeping K\ fixed, K 2 can be chosen in such a 
way to make the second term as small as we want too. At last, keeping K 2 fixed, one can 
choose K 3 large enough to make the third term as small as we want (as can be computed 
since the limit is given by the K 3 correlation function of the Airy kernel). □ 



To complete the proof of Hypothesis |4.1[ we therefore need to show that 



Lemma 5.6. Assume that the entries of X n satisfy Hypothesis |5.ffi Then, for any 5 > 0, 
any finite integer number p, 



lim P 

n— >oo 



n 

- T 

n A" + 2 

3=p+l J 



1 



> 5 







Proof. Notice that by (HI EH] we know that the p smallest eigenvalues of X n converge 
in law towards the Tracy-Widom law, so that 



lim lim P ( min |A™ + 21 < en 2/3 



n— >co \ l<j<p 

Thus, for any finite p, with large probability, 

1 



0. 



n<-^ A™ + 2 

3=2 J 



< 



pe n a 



and therefore it is enough to prove the lemma for any particular p. As in the previous 

2 

proof, we choose p large enough so that A" > —2 + with probability greater than 
1 — 8{p) with 5(p) going to zero as p goes to infinity. We shall prove that with high 
probability 



j [7"] j 

lim lim — > 

74,0 rwoo n ^ A™ + 2 

j=p J 



< 0. 



(43) 



This is enough to prove the statement as for any 7 > 0, 2 + \? , converges to 5(7) > so 
that fix([8(j), 2]) = 1—7, see [331 Theorem 1.3], 



lim 



n 

- E 



[nrf\ 



-5(7) 



2 + X 



which converges as 7 goes to zero to J (2 + x) 1 dfix(x) = 1 (by e.g. (37)). To prove (43) 



we choose p G (2/3, a/2/3) and write, on the event A™ + 2 > Ap + 2 > n 3 > n p for 
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3 > P, 



hn] 



n 4^ A? + 2 



b/n] i K 



Kk<K 



J=2 



n(A? + 2) 



For the first term, we use Sinai-Soshnikov bound, which under the weakest hypothesis are 
given in [39, Theorem 2.1]. It implies that with probability going to one with M going to 
infinity, for s n = o(n 2//3 ) going to infinity, 



E 

i=i 



X? 



< M 



n 

~ ■ 

o n 



i fc+i 

This implies, by Tchebychev's inequality and taking s n = n +p that 



n r , n 









2 



> 1 - n" 



< (1 - n" 



n fc+i , 



E 

i=i 



K 

2 



3,k+l 



Consequently we deduce that 

A n <eM ^ 



n fc —I fc+1 



rr n 2 



l<fe<X 

which goes to zero as p > 2/3. For the second term B n , note that by [4*21 Theorem 1.10], 
for any e > small enough, 

\N n [n~ e i, n~ e {£ + 1)] - np x {[-2 + n~% -2 + n^(£ + 1)]) | < n 1 - 5 ^ 

with 5(e) = Hence, since [ix{[— 2 + n~ e £, —2 + n~ e (£ +1)]) ~ V^, we deduce for 
e small enough that for all I > 1, 

A4K e £,n~ £ (£+ 1)] < 2n l -'^Vl. 

This allows to bound B n by 



5„<2V ^-n - ^ v 7 ^ < 2 / — = 2^ 
£=1 * Jo V x 



which goes to zero as n goes to infinity and then 7 goes to zero. 



□ 



5.2. Coulomb Gases. We can also consider random matrices X n which law is invariant 
under the action of the unitary or the orthogonal group and with eigenvalues with law 
given by 

n 

dP n (X 1 ,...,X n ) = ^|A(A)|^-^^=x^)TTdA, (44) 

i=l 

with a polynomial function V of even degree and positive leading coefficient and /3 = 1, 2 
or 4. We assume moreover that V is such that the limiting spectral measure \xy of (X n ) 
is connected and compact and that its smallest and largest eigenvalues converge to the 
boundaries of the support. This set of hypotheses is often referred to as the "one-cut 
assumption" . It holds in particular if V is strictly convex and this includes the classical 
Gaussian ensembles GOE and GUE (with V(x) = x 2 /4 and (3 = 1,2). 
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Proposition 5.7. Under the above hypothesis on V, the extreme eigenvalues of X n con- 
verge to the boundary of the support. The convergence of the extreme eigenvalues of X n is 



given by Theorem 2.1. These eigenvalues have Gaussian fluctuations as stated in Theorem 



3.2 if they deviate away from the bulk. 



Suppose moreover that Assumption ^. 2\ holds. 

If the perturbation is of rank one and is strong enough so that the largest eigenvalues 
deviates from the bulk, for all k > 2, the rescaled kth largest eigenvalue ns(A"_ fc+1 — by) 
converges weakly towards the k — 1-th Tracy Widom law. If the perturbation is of rank 
one and is weak enough, for all k > 1, the rescaled kth largest eigenvalue ns(A^_ fc+1 — by) 
converges weakly towards the k-th Tracy Widom law. 

If the perturbation is of rank more than one, the extreme eigenvalues of X n sticking to 
the bulk are at distance less than n~ 1+e for any e > from the eigenvalues of X n . In the 
i.i.d model, Theorem \4-5\ prescribes exactly in the neighborhood of which eigenvalues of 
X n each of them lie. 

Proof. As explained above, it suffices to verify that the hypotheses hold in probability 
for {X n ) n > x . 

Note that the convergence of the spectral measure, of the edges and the fluctuations 
of the extreme eigenvalues were obtained in [47J. The fact that y/n{G^ n {z) — G sc (z)) 



converges in probability to zero is a consequence of [28] so that Hypothesis 3.1 holds 



a > 1/3 and any integer p. We first show that 



We next check Hypothesis 4.1[p,a] for the matrix model P n . We shall prove it for any 

-V'(a v ) . 



lim E 

n— >oo 



(45) 



Indeed, the joint distribution of (A™, . . . , A") is 

1 - 

e-nE^vW Y[ (Aj — Aj) /3 lA n dAi • • • dA n , 

n l<i<j<n 

with (3 = 1, 2 or 4, Z@ is the normalising constant and A n = {Ai < • • • < A n }. 
Therefore, 



E 



J i\ n P i<i<j< n 

^1 OX ll A: ) ft (A l -A,)^dA 1 ...dA n 



-n/3E [V'(X n p )] , 



by integration by parts. Equation (|45|) follows, since A^ converges almost surely to ay 



(and concentration inequalities insures V'(Xp) is uniformly integrable). But, for any e > 0, 
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with, by convergence of the spectral measure and of A™, the right hand side converging 
to — G^ x (— ay — e) which converges as e decreases to zero to — G llx {— ay) = —V'(ay). 
Hence, ^ J2i^ P a"— A n * s bounded below by —V'(ay) with large probability for large n, and 
converges in expectation to —V'(ay), and therefore converges in probability to —V'(ay). 
Moreover, by jH] (see [3S] in the Gaussian case), the joint law of 

(n 2 /\Xl - a v ),n 2 '\\ n 2 - a v ), . . . ,n 2 / 3 (A£ - a v )) 

converges weakly towards a probability measure which is absolutely continuous with re- 
spect to Lebesgue measure. As a consequence, we also deduce from the first point that 
rT x ^2 i<m (A™ — A™) -1 vanishes as n goes to infinity in probability for m n <C n 1//3 and 



therefore (|45j) proves the lacking point of Hypothesis 4.1 

For the two other points, observe that [37] implies that for any e > 0, P(|AJj — A"| < 

2 2 2 

n~*~ e ) — Y 0. On the event {\X% — A™| > 7?,~3~ e }, we have |A™ — A™| > n~s~ e for all 
i £ [2, n — 1], so that 

— > <n~* +e -y 



n2 U ( A " - A ") 2 " n U A " - A " 

JL V 1 < n- 1+3e - V 1 

n 4 (K ~ A™) 4 " n ^ A" - X'i 



so that by (45) and Markov's inequality, Hypothesis 4.1 holds in probability for any 



77 < 1/3, 774 < 1 and a > 1/3. □ 

5.3. Wishart matrices. Let G n be an n x m real (or complex) matrix with i.i.d. centred 
entries with law /1 such that f zdfi(z) = 0, f \z\ 2 dfi(z) = 1 and f \z\ 4 dfi(z) < 00. Let 
X n = G n G*Jm. 

Proposition 5.8. Let n,m tend to infinity in such a way that n/m — > c G (0,1). The 



limits of the extreme eigenvalues of X n are given by Theorem 2. 1 and the fluctuations 



of those which limits are out of [a, b] are given by Theorem 3.2, where the parameters 



a, b, p 9 , c a are given by the following formulas: a = (1 — \/c) 2 , b = (1 + a/c) s 




e-c 



if \9 — c\ > y/c, 

if \9 — c\ < y/c and 9 > 0, 

if \9 — c\ < y/c and 9 < 0, 



and 



a<2 fl — (a-c) 2 ^) ^ n the i.i.d. model, 



1 1 — 7 — ) in the orthonormalised model. 



Assume now that the law of the entries satisfy Hypothesis 5.2 If the perturbation has 
rank one, we have the following precise description of the fluctuations of the extreme 
eigenvalues of X n : 
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• If 6 > c + a/c (Vesp. 9 < c — yfc), for all p > 2, n 2 / 3 (A™„ p+1 — 2a) (resp. 
n 2//3 (A™ — 2a) ) converges in law to the p — 1th Tracy Widom law. 

• //0<6'<c+ v /c (resp. c — a/c < 6 < 0), for all p > 1, n 2 / 3 (A™_ p+1 — 2a) (Vesp. 
n 2//3 (A™ — 2a) ) converges in law to the pth Tracy Widom law. 

If the perturbation has rank more than one and for alii, 9i {c+ a/c, c— a/c}, t/ie extreme 
eigenvalues of X n are at distance less than n~ 1+e for any e > to the extreme eigenvalues 
of X n , which have Tracy- Widom fluctuations. 

Before getting into the proof, let us make a remark. The Proposition above generalizes 
some results first appeared in [91 [19] . In these papers, the authors consider models with 
multiplicative perturbations (in the sense that the population covariance X matrix is 
assumed to be a perturbation of the identity). Here, we consider additive perturbations 
but the two models are in fact similar, since a Wishart matrix can be written as a sum of 
rank one matrices YlT=i <J i^i^i i with <jj the eigenvalues of £ and Y,i n-dimensional vectors 
with i.i.d. entries. So, adding our perturbation Y^i=i @iUiU* boils down to change m into 
m + r (the limit of m/n is not changed) and to extend S with some new eigenvalues 
#i, . . . , 6 r . 

Proof. Again, it suffices to verify that the hypotheses hold in probability for (X n ) n >i. 

It is known, [32], that the spectral measure of X n converges to the so-called Marcenko- 
Pastur distribution 

dp:x{x) ■= w—V( b -x){x- ajl [ab] (x)dx, 

where a = (1 — a/c) 2 and b = (1 + a/c) 2 . It is known, [U Th. 5.11], that the extreme 
eigenvalues converge to the bounds of this support. The formula 



n ( \ z + c ~ l ~ s S n - a ) V( z - c - l) 2 - 4c 

GW (z) = i^f [zeR\ [a, b\ ) 

allows to compute p e and c a . Moreover, by [3J Th. 1.1] or [U Th. 9.10], we also know 
that a central limit theorem holds for the linear statistics of Wishart matrices, giving 



Hypothesis 3.1 as in the Wigner case. 



For Hypothesis 4.1, the proof is similar to the Wigner case. The convergence to the 
Tracy- Widom law of the non-perturbed matrix is due to S. Peche [37] (see [33J and [2D J for 
the Gaussian case). The approximation of the eigenvalues by the quantiles of the limiting 
law can be found in [T7J Theorem 9.1] whereas the absolute continuity property needed 



to prove Lemma 5.5 is derived in [T71 Lemma 8.1]. This allows to prove Hypothesis 4.1 



in this setting as in the Wigner case, we omit the details. □ 

5.4. Non-white ensembles. In the case of non-white matrices, we can only study the 
fluctuations away from the bulk (since we do not have the appropriate information about 



the top eigenvalues to prove Hypothesis 4.1). We illustrate this generalisation in a few 



cases, but it is rather clear that Theorem 3.2 applies in a much wider generality. 



5.4.1. Non-white Wishart matrices. The first statement of Proposition 5.8 can be gener- 

1 1/2 1/2 1 

alised to matrices X n of the type X n = —T n G n G* n T n or —G n T n G* n1 where G n is annxm 
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real (or complex) matrix with i.i.d. centred entries with law \i such that J zdfj,(z) = 0, 
J \z\ 2 dfi(z) = 1 and J \z\ 4 dn(z) < oo and T n is a positive non random Hermitian n x n 
matrix with bounded operator norm, with a converging empirical spectral law and with 
no eigenvalues outside any neighborhood of the support of the limiting measure for suffi- 
ciently large n. Indeed, in this case, everything, in the proof, stays true (use [21 Th.1.1] 
and jH Th. 5.11]). However, when the limiting empirical distribution of T n is not a Dirac 
mass, the computation of the p^'s and the c a 's is not easy 

5.4.2. Non-white Wigner matrices. There are less results in the literature about the cen- 
tral limit theorem for band matrices (with centring with respect to the limit) and the con- 
vergence of the spectrum. We therefore concentrate on a special case, namely a Hermitian 
matrix X n with independent Gaussian centred entries so that i^flX^! 2 ] = n~ 1 <j{i/n,j/n) 
with a stepwise constant function 

k 

In [3T], matrices of the form S n = Ylj=i a i ® w hh some independent matrices Xj 
from the GUE and self-adjoint matrices aj were studied. Taking a,j = (e Py e + £e,p)&p,e or 
i( e p,e ~ £l,p)&p,t with e Py e the matrix with null entries except at (p,£) and 1 < p < £ < k, 
we find that X n = S n . Then it was proved [311 (3.8)] that there exists a, e, 7 > so that 
for z with imaginary part greater than n -7 for some 7 > 0, 



E 



n 



-Tr(z-X n ) 



G{z) 



a n -l-e 



(46) 



which entails the convergence of the spectrum of X n towards the support of the limiting 
measure [31, Proposition 11] with exponential speed by [3T1 Proof of Lemma 14]. Thus 
X n satisfies Hypothesis 1.1 Hypothesis |3.1 can be checked by modifying slightly the proof 



of (46) which is based on an integration by parts to be able to take z on the real line 



but away from the limiting support. Indeed, as in [23], Section 3.3], we can add a smooth 
cut-off function in the expectation which vanishes outside of the event A n that X n has 
all its eigenvalues within a small neighborhood of the limiting support. This additional 
cut-off will only give a small error in the integration by parts due to the previous point. 



Then, (46), but with an expectation restricted to this event, is proved exactly in the same 
way, except that Qz can be replaced by the distance of z to the neighborhood of the 
limiting support where the eigenvalues of X n lives. Finally, concentration inequalities, in 
the local version [221 Lemma 5.9 and Part II], insure that on A n , 



1 



-Tr(z - X ri 



E 



n 



l An -Ti(z - X n 
n 



is at most of order n 1+6 with overwhelming probability. This completes the proof of 



Hypothesis 3.1 



5.5. Some models for which our hypothesis are not satisfied. 
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We gather hereafter a few remarks about some models for which the hypothesis we 
made on X n are not satisfied. For sake of simplicity, we present hereafter only the case 
of i.i.d. perturbations (1). 

5.5.1. I.i.d. eigenvalues with compact support. We assume that X n is diagonal with i.i.d. 
entries which law p is compactly supported. As in the core of the paper, we denote by 
a (resp. b) the left (resp. right) edge of the support of p. We also denote by its 
cumulative distribution function and assume that there is k > such that for all c > 0, 

hm ^- f = c K (47) 

*->o+ 1 - F^b-x) y 1 

In this situation, it is easy to check that Hypothesis |1.1| holds in probability with 



Px = P- But Hypothesis |3.1| is not satisfied. Indeed, by classical CLT, we have, for 
Pa £ [a, 6], 

Wl = ^a{G,M ~ G„(p«)) 
converges in law, as n goes to infinity to a Gaussian variable W a with variance —G' (p a ) — 
G>(p a ) 2 . Moreover, 

E[W a W a ,} = [ - -i -dfi(X) - G^ Pa )G,{ Pa ,). 

J [Pa — A){Pa> — 



Nevertheless, Theorem [372] holds for this model. Indeed, the whole proof of this theorem 
goes through in this context, except the proof of Lemma [3. 5[ where we have to make the 
following decomposition M™ t (i, x) = M™f(i, x)+M?f(i, x) + M"f (z, x) with the difference 
that this time M™f does not go to zero but converges towards W ai . Hence, the eigenvalues 
fluctuate according to the distribution of the eigenvalues of (cjMj + W aj Ikj)i<j< q , with Cj 



and Mj as in the statement of Theorem 3.2 and Ik denotes the kj x kj identity matrix. 



Let us now consider the fluctuations near the bulk. We first detail the fluctuations of the 
extreme eigenvalues of X n . According to [26], the fluctuations of the largest eigenvalues 
of X n are determined by the parameter k defined in (47), that is, if v n = F^ip — 1/n), 



then the law of converges weakly to the law with density proportional to e~ x on M . 
Otherwise stated, the fluctuations of A" are of order rT x l K with asymptotic distribution 

rat if k < 1, then 6 = 0. 

p, a] to hold, we need a > \ — \ 



4.1 



the Gumbel distribution of type 2. One can check t 
One can show that, for any fixed p, for Hypothesis 
and we then obtain that the distance of the extreme eigenvalues of the deformed matrix 
is at distance less that n~ 1+a for any a' > a. Therefore if k > 4/3, this theorem allows 
us to deduce that the fluctuations of the extreme eigenvalues of the deformed matrix are 
the same as those of the non-deformed matrix. 

5.5.2. Coulomb gases with non-convex potentials. In [32], Pastur showed that for a Coulomb 



gas law (44) with a potential V so that the equilibrium measure has a disconnected sup- 
port, the central limit theorem does not hold in the sense that the variance may have 
different limits according to subsequences (see [351 (3-4)]. Moreover the asymptotics of 
^/n(TY(X n ) — p(x)) can be computed sometimes and do not lead to a Gaussian limit. 
We might expect then that also \fn{G IJin {x) — G^(x)) converges to a non-Gaussian limit, 
which would then result with non-Gaussian fluctuations for the eigenvalues outside of the 
bulk. 
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6. Appendix 

6.1. Determinant formula. We here state formula (|3]), which can be deduced from the 
well known formula det ( ^ ^ J = det(D) det(A — BD~ l C). 

Lemma 6.1. Let z G C\{A", . . . , A™} and 9i, . . . , 9 r ^ 0. Set D = diag(6 l 1 , . . . , 6 r ) and 
let V be any n x r matrix. Then 

det (z-X n - VDV*) = det(z - X n ) det(D) det (IT 1 - V*(z - X n )~ V) 

6.2. Concentration estimates. 



Proposition 6.2. Under Assumption t/iere exists a constant c > so that for any 
matrix A := (djk)i<j t k<n with complex entries, for any 5 > 0, for any g = (gx, . . . ,g n ) T 
with i.i.d. entries (<?i)i<i< n with law u, 

P (\(g, Ag) - E[(g, Ag)} \ > 6) < Ae~ c ^^ } 
if C 2 = Tt(AA*) and if g is an independent copy of g, for any 5, k > 0, 

P (\{g,Ag)\ > 5\Jtt(AA*) + Ky/Ti((AA*j^J < 4e~ cfi + 4 e - cmin {^ 2 }. 

Proof. The first point is due to Hanson- Wright Theorem [21], see also [131 Proposition 
4.5]. For the second, we use concentration inequalities, see e.g. [H Lemma 2.3.3], based 
on the remark that for any fixed g, g — )■ (g, Ag) is Lipschitz with constant a/ (g, AA*g) 
and therefore, conditionally to g, for any 5 > 0, 

P (\{g,Ag)\ > 8^/(9, AA*~g)) < Ae~ c&2 

On the other hand, the previous estimate shows that 

P (\(g,AA*g) - Tr(A4*)| > k^Tt(AA*) 2 ^ < 4 e ~^{^ 2 } . 

As a consequence, we deduce the second point of the proposition. □ 

Let G n = [g™ ■ ■ ■ g™\ be an n x r matrix which columns g™, . . . , g™, are independent 
copies of an n x 1 matrix with i.i.d. entries with law v and define 

V^ = -(g?,g]), l<i,j<r, 

and, for j < i - 1, if det^]^ ^ 0, 

_ detbff]^ • _ f V%, if I * j, 

On det[V^i]]^ =1 = 0, we give to W£j an arbitrary value, say one. Putting W£ = 1 and 
W?j = for j > i + 1, it is a standard linear algebra exercise to check that the column 
vectors 



^ W^g] = zth column of G"(H/ n ) T 

are orthogonal in C n . Let us introduce, for M an r x r matrix, ||iW = sup 1< j J<r 
We next prove 
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Proposition 6.3. For any 7 > 0, there exists finite positive constants c,C (depending on 
r) so that for Z n = V n or W n , 



+ e" 



Moreover, with ||i>||| — Y^=i \ v i\ 2 > f or an V 7 e (0? V^(2 r — e) for some e > 0, 
1 



P I max 

Ki<r 



> n-27 < C 



-4-i C 2-V +4e - 



Proof. We first consider the case Z n = V n . The maximum of \V£j — 5y \ is controlled by 
the previous proposition with A = n~ x l , and the result follows from TrAA* = rC 1 and 
Tr((AA*) 2 ) = n~ 3 , and choosing 5 = 7/v^2, k = Vn. The result for W n follows as on 
HV^-JHoo <7n"5 < 1 

ldet^Lx-11^2^-5, 

whereas 

|det[7S'']£i=il^ 2P 'm-*. 
For the last point, we just notice that since ^|| Y7j=i ^'ijQj'Wl = (ZVZ*) iti , we have 



max 

Ki<r 



'iiE z .wii2-i 



n 



3=1 



<C(r) max ||Z n |P max \\Z n - JlL 

Z"=P or H 7 " Z n =V n or W n 



for a finite constant C(r) which only depends on r. Thus the result follows from the 
previous point. □ 



6.3. Central Limit Theorem for quadratic forms. 



Theorem 6.4. Let us fix r > 1 and let, for each n, A n (s,t) (1 < s,t < r) be a family 
of n x n real (resp. complex) matrices such that for all s,t, A n (t,s) = A n (s,t)* and such 
that for all s,t — 1, . . . ,r, 



in the i.i.d. model, 

-Tr[A n ( S ,t)A n (s,tT} -^al, ±J2\A n (s,s) iti \ 



(48) 



1=1 



in the orthonormalised model, 



-Tr[\A n (s,t)--TrA n (s : t)\ 2 } — ► a% t 

n ?i-5>oo n 



n 



1 n 



i=l 



A n {s,s\ i - -Tr A n (s,t) 



(49) 



for some finite numbers o- s j,u s (in the case where k±{v) = 0, the part of the hypothesis 
related to u s can be removed). For each n, let us define the r x r random matrix 



s/n [ {u:,A n (s,t)u?) - l s=t -Tr(A n (s,s)) 

71/ 



s,t=l 



EXTREME EIGENVALUES OF DEFORMED RANDOM MATRICES 



37 



Then the distribution of G n converges weakly to the distribution of a real symmetric (resp. 
Hermitian) random matrix G = [g s ,t] r s t=i suc -h that the random variables 

{g s ,t ;i<s<t<r} 

(resp. {g s , s ; 1 < s < r} U {^(g s , t ) ; 1 < s < t < r} U {%(g s , t ) ;l<s<t<r}) 

are independent and for all s, g S;S ~ Af(0,2a^ s + /t 4 (z/)u; s ) (resp. g StS ~ A/"(0,cr 2 s + 
Ki(v)uj s )) and for all s ^ t, g s>t ~ jV(0, a 2 s>t ) (resp. $t(g S}t ), $S{g,,t) ~ AT(0, <7 2 t /2)). 

Remark 6.5. iVo£e i/iai z/ £/ie matrices A n (s,t) depend on a real parameter x in such a 
way that for all s, t, for all x, x' G R, 



0. 



7? 



then it follows directly from Theorem 6J^_ and from a second moment computation that 
each finite dimensional marginal of the process 

fn-[(u n s1 A n {s,t){xsM) - l s=t -Tr(A n ( S , S )(x s , s )) 



n 



l<s,i<r, x s , t eR, x Stt =x tjS 

where there is no 



converges weakly to the law of a limit process [g s ,t]i<s,t<r , x s t &R, x s t =x t 
dependence in the variables x s>t (1 < s,t < r). 

Proof. • Let us first consider the model where the (y/nu™ )i< s < r are i.i.d. vectors with i.i.d. 
entries with law v satisfying Assumption 1.2 Note that for all s,t = 1, . . . , r, by (48), 
the sequence ^ Y^j=i A n (s, is bounded. Hence up to the extraction of a subsequence, 
one can suppose that it converges to a limit r St t G C. Since the conclusion of the theorem 
does not depend on the numbers r Sj t and the weak convergence is metrisable, one can 
ignore the fact that these convergences are only along a subsequence. In the case where 
Kn{v) = 0, we can in the same way add the part of the hypothesis related to co s . 

We have to prove that for any real symmetric (resp. Hermitian) matrix B := [& s ,t]s£=i, 
the distribution of Tr(BG n ) converges weakly to the distribution of Tr(BG). Note that 

Tr(BG n ) = -^(U*C n U n -TrC n ), 



n 



where C n is the rn x rn matrix and U n is the rnxl random vector defined by 



b ltr A n {l,r)' 
b, 



u„ 



b r ,iA n (r, 1) ■■■ b r>r A n (r,r)_ 

In the real (resp. complex) case, let us now apply Theorem 7.1 of [7] in the case K — 1. 
It follows that the distribution of 



Tr(BG n ) = J2KsG n ,s,s+ ^(KtMGn,s,t) + 2Q(b S}t MG 

l<s<t<r 



n,s,t , 



converges weakly to a centred real Gaussian law with variance 



l<s<t<r 



It completes the proof in the i.i.d. model. 



in the real case, 

in the complex case. 
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• In the orthonormalised model, we can write u" = -n^s — \, m M , W^Oj, where the 

s \\h,i=i W^gi\\2 ^-O-i s 3VJi 

matrix W n is the one introduced in this section. It follows that, with 



B n (s,t) = A n (s,t) - ^Tr(A n (s,t)), 



by orthonormalization of the u"'s 



n( «,A"(,,tK)-— Tr(A"( S ,t)) 



n 



n(ii;,B»(MK) 



t 
n 



But, by the previous result, if i ^ j, 

(gi,B(s,t)gj) 



converges in distribution to a Gaussian law, whereas if % = j, 

(9i,B(s,t)gi) 



n 



1 Ty(A(s t)) 

= {{ 9i , A(s, t) 9i ) - E[(g h A(s, t) gi )\) + 1 y " {{g h gi ) - E[{ 9i , 9i )]) 



where both terms converge to a Gaussian. Thus this term is also bounded as n goes to 
infinity. 



Hence, by Proposition 6.3, we may and shall replace W n by the identity (since the error 



term would be of order at most n 2 +e ), which yields 

v^«,ir( s ,tK> ~ V^ l {9s,B( Sl t)g t ) 

so that we are back to the previous setting with B instead of A. □ 

Acknowledgments: We are very grateful to B. Schlein for communicating us Lemma 
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